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Abstract: Hypermetabolism is a significant sequela to severe trauma such as burns, as well 
as critical illnesses such as cancer. It persists in parallel to, or beyond, the original 
pathology for many months as an often-fatal comorbidity. Currently, diagnosis is based 
solely on clinical observations of increased energy expenditure, severe muscle wasting and 
progressive organ dysfunction. In order to identify the minimum number of necessary 
variables, and to develop a rat model of burn injury-induced hypermetabolism, we utilized 
data mining approaches to identify the metabolic variables that strongly correlate to the 
severity of injury. A clustering-based algorithm was introduced into a regression model of 
the extent of burn injury. As a result, a neural network model which employs VLDL and 
acetoacetate levels was demonstrated to predict the extent of burn injury with 88% 
accuracy in the rat model. The physiological importance of the identified variables in the 
context of hypermetabolism, and necessary steps in extension of this preliminary model to 
a clinically utilizable index of severity of burn injury are outlined. 

Keywords: hypermetabolism; metabonomics; flux analysis; cluster analysis; regression; 
variable selection; burn injury 



Metabolites 2012, 2 



459 



1. Introduction 

Hypermetabolism is a significant consequence of severe trauma such as burns [1], as well as critical 
illnesses such as cancer [2]. It persists in parallel to, or beyond, the original pathology for many 
months as an often-fatal comorbidity. Prolonged hypermetabolism is characterized by increased resting 
energy expenditure and severe muscle wasting due to a negative nitrogen balance. The underlying 
mechanisms that control the onset and resolution of the hypermetabolic response are unknown. Therefore, 
current treatments are directed at symptoms which are metabolic, endocrine, and immune in nature. 

Increasing nutritional energy delivery and protein intake only partially alleviate the loss of lean 
body mass [3,4]. Experimental approaches to overcome the deleterious effects of hypermetabolism 
have been used with varying success, including glutamine and arginine supplementation [5,6]; 
combinatorial nutritional therapies using a diet high in vitamins, protein, amino acids, and co-3 fatty 
acids [7]; peroxisome proliferator activated receptor-a agonists to improve fat oxidation and 
mitochondrial activity [8]; as well as antioxidant and anti-inflammatory agents [9]. Modulation of 
insulin action [10], direct insulin therapy [11,12], administration of other anabolic agents [13], as well 
as P-blockers [14] have also produced significant improvements, but are inherently impractical in the 
long-term, and in some cases produce undesirable, potentially fatal, side effects [15,16]. 

Rational design and optimization of nutritional therapies can be achieved by targeting the 
interconnected metabolic network and regulatory pathways impacted by hypermetabolism [17-21], 
thereby having a broad impact on the metabolome, and subsequently altering the physiome at the 
genomic and proteomic levels. Hence an understanding of the connections among individual 
metabolites and the overall injury physiome is necessary to rationally design metabolic interventions. 
While the advances in metabolomics and metabonomics are starting to present an ever-increasing 
amount of data, the capability to quantitatively and accurately predict the metabolic effects of 
nutritional supplements in order to design the kind of combination therapies that could treat 
hypermetabolism and other metabolic conditions remains absent. 

There are a variety of techniques to extract knowledge from "omics" data. These include pattern 
identification methods such as clustering and principal component analysis (PCA) typically applied to 
time-series mRNA microarray data to identify key trendlines [22]; network analysis to identify 
correlations between genes/metabolites [17]; and marker discovery [23], which utilizes a combination 
of techniques above to correlate disease with a specific measured variable. Metabonomics in particular 
focuses on the metabolic analysis of the consequences of a perturbation, such as disease and 
medication. While metabonomics of toxicity [24] and pharmaceuticals [25] is an emerging field, there are 
few studies performing rigorous metabolic analyses as a function of extent of injury or disease [23,26]. 

While there has been some effort to develop analytical mathematical models for hypermetabolism 
[27,28], predictive models remain absent. There are however efforts to develop such approaches, for 
instance Flux Balance Analysis (FBA), to identify the metabolic causes of fat accumulation in 
hepatocytes [29,30]; similar approaches could potentially predict metabolic responses to potential 
interventions, such as amino acid supplements and hence rationally develop combination therapies in 
hypermetabolism. However, in vivo hypermetabolism remains a mostly clinical observation, and 
quantitative measures of hypermetabolism to define the extent of burn injury are necessary for use of 
optimization-based FBA models. 
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The objective of this study was to identify sensitive indicators that correlate with the severity of 
injury and that can be measured from blood samples. For this purpose, we used rat models of 
cutaneous burn injury of increasing size, and analyzed metabonomic data on postburn day 4. We 
identified VLDL and acetoacetate in the circulation as being particularly sensitive to the extent of burn 
injury, and thus could serve as a quantitative measure of the grade of hypermetabolism. 

2. Experimental Methods 

Briefly, male Sprague-Dawley rats (Charles River Labs) weighing between 270 g and 300 g at time 
of burn, were subjected to a third degree cutaneous burn injury (i.e., depth of injury spans the entire 
thickness of the skin) covering 20% of the Total Body Surface Area (TBS A) (dorsal burn only) or 40% 
TBSA (burn on dorsum and abdomen) by contacting the skin with water at 100 °C. Sham-treated 
animals were handled identically to the burn groups, except that room temperature water was used. 
After injury, animals were immediately resuscitated with an intraperitoneal saline injection (20 mL/kg 
per rat) and allowed to recover in individual cages. Animals were weighed daily and food consumption 
was monitored. On the third day following the injury, all rats were fasted overnight in preparation for 
the blood samples to be taken on Day 4 as described in detail elsewhere [21]. Briefly, on day 4, each 
rat was anesthetized and blood flow through each of the major blood vessels entering the liver 
(the portal vein, PV, and hepatic artery, HA, USA) were measured using a perivascular ultrasonic 
flow-probe (Transonic Systems, Ithaca, NY, USA). The sum of flow rates into the liver was assumed 
to equal the flow rate out of the hepatic veins into the suprahepatic vena cava (SHVC). Following flow 
rate measurements, blood samples from the hepatic veins and PV were taken, followed by arterial 
blood. Blood samples were analyzed for blood gases and pH using a Rapidlab Blood Gas Analyzer 
865 (Bayer). Standard reagent kits were used to determine plasma glucose (Stanbio No. 1075-825), 
urea (Stanbio No. 0580), and lactate (Trinity Biotech USA No. 735-10). The ketone bodies 
acetoacetate and [3-hydroxybutyrate were measured enzymatically by following the appearance or 
disappearance of NADH upon the addition of [3-fiydroxybutyrate dehydrogenase (Sigma), respectively. 
Nineteen of the common amino acids (except tryptophan) plus ornithine and ammonia were measured 
using a Waters HPLC apparatus (Waters Co. Milford, MA, USA). ELISA techniques were used to 
detect albumin (Sigma) and insulin (Crystal Chem Inc No. INSKR020). Alkaline phosphatase (ALP), 
alanine transaminase (ALT), aspartate aminotransferase (AST), total bilirubin, blood urea nitrogen, 
and creatinine, HDL, LDL, VLDL, cholesterol and triacylglycerols were measured using a Piccolo 
Comprehensive Metabolic Panel (Abaxis, Inc., Union City, CA, USA). Finally, the liver was excised 
and weighed. Fluxes across the liver were subsequently determined as the difference between in- and 
effluxes, calculated per vessel as the product of metabolic concentration and flow rate normalized to 
the weight of the liver. Note that early injury markers, such as cytokines and tissue damage, have 
shown a clear linear trend with increasing burn severity [31,32]; however, we sought variables that 
were less likely to have a transient response or a response subsequently complicated by sepsis so we 
only employed metabolic markers still visible 4 days after injury for consideration as indicators of 
burn injury. 
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3. Numerical Methods 

3.1. Theoretical Aspects 

The primary goal of this study was to identify key metabolites that are strongly correlated to the 
hypermetabolism in the form of a predictive model of the degree of burn injury. However, since 
metabolic data suffer from the existence of a high degree of correlation [17], a straightforward 
statistical analysis (e.g., ANOVA), is very limited. Many variables will be identified as correlated to 
injury because of their own interrelationships. Therefore, a critical issue is to identify the key variables 
involved, which may be defined as the minimum number of variables that can explain the metabolic 
response to injury. For the purposes of this work, this problem can equivalently be stated as the 
identification of the minimum number of variables necessary to construct an accurate model that can 
predict the degree of injury from metabolic data. 

It should be noted that in the absence of a mechanistic description of the effects of injury, empirical 
mathematical models are necessary, which further increases the problem of complexity as the type of 
model to be used becomes another variable. Therefore, the problem of constructing an "index of burn 
injury severity" is a multiobjective problem where the task is to simultaneously: (i) select and train the 
best mathematical model; (ii) maximize model accuracy by selecting the variables to use; (iii) minimize 
the number of variables in the model. To achieve this goal, we designed a novel algorithm to identify the 
metabolites that are most indicative of injury grade, as outlined in Figure 1 and discussed in detail below. 

Figure 1. The burn-injury indicator identification algorithm. 
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identify patterns 
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the best classifier with N variables 

Clustering of Dose-Response Patterns. The first step in the algorithm is clustering in order to group the 
metabolites responding similarly to increased injury, which serves multiple purposes. Clustering serves 
as the first step in reducing the problem dimensionality, as all metabolites collected in one group can 
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be considered to be regulated by the same mechanisms, and ideally can be considered as a single 
variable (or, a single metabolite can adequately represent all others in the same cluster). Therefore, the 
number of clusters necessary to represent the data presents an adequate first-guess for the ideal number 
of variables to capture the entire metabolic response. In addition, since we have obtained 
dose-response data rather than a simple injury vs. non-injury comparison, clustering also identifies the 
major patterns observed in response to increasing levels of burn injury, which enables incorporating a 
physiological interpretation prior to identification of the burn-grade indicators. 

A critical issue in clustering is the determination of the number of clusters, either defined directly as 
a parameter (as in k-means clustering which is employed here) or indirectly (as in hierarchical 
clustering). We utilized a tandem approach to determine the optimum number of clusters, combining 
analysis of explained variance (here measured as sums of point-to-centroid distances of all clusters) 
and the separation of individual clusters (assessed through silhouette analysis) [33]. It should be noted 
that the tandem use is necessary: in clustering, by definition, as the number of clusters increases, the 
error of clustering is reduced, until each instance {i.e., patient) itself is a cluster and the error is zero, 
which obviously is of little value. An alternative criterion to identify a good number of clusters that is 
commonly employed is the marginal gain of adding a cluster, which often displays "elbows" where the 
marginal value decreases steeply. However, often there are multiple elbows and so multiple choices of 
good numbers of clusters. To differentiate between such potential solutions, we employed the 
silhouette analysis, which is a measure of separation for clusters. Briefly, the silhouette value is a 
measure of how close each point in one cluster is to points in the neighboring clusters. This measure 
ranges from +1, indicating points that are very distant from neighboring clusters, through 0, indicating 
points that are not distinctly in one cluster or another, to — 1, indicating points that are probably 
assigned to the wrong cluster [33]. Combining the marginal value and silhouette methods, it is possible 
to assess marginal return and mean silhouette values as a function of number of clusters, and identify 
the number of clusters that simultaneously have a local maxima in their silhouette value and display a 
decrease in cluster error. This is in essence a visual analysis of the two charts to identify the minimum 
number of clusters 

Identification of Independent Patterns. Identified patterns reveal the key trends in the dose response. In 
mRNA analysis, the co-regulated gene expression reveals information regarding genetic control 
motifs [34]. The situation is more complex in metabonomics, as correlations exist due to a variety of 
reasons such as simple stoichiometric dependencies [17]. Accordingly, it is necessary to extract the 
actual number of independent patterns in order to filter out these metabolic correlations that are not of 
primary interest for this study. 

To achieve this purpose, a Singular Value Decomposition (SVD) was performed on the cluster 
centroids: Briefly, each cluster was expressed as a vector (the cluster median for each level of injury 
was an entry in the vector), and the vectors were augmented to form a cluster centroid matrix. SVD 
was performed to identify the number of independent patterns, which can be used to explain the 
remaining patterns. This is in effect a rank analysis, but direct analysis of singular values enables 
identification of marginally non-zero singular values, which are likely artifacts of measurement error 
and/or statistical clustering rather than major patterns in the injury response. The number of distinct 
singular values that are also non zero (>10 , a heuristically chosen limit) was chosen as the number of 
independent clusters. 
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The product of this tandem approach is: (i) identified patterns in the injury response; and (ii) the 
minimum number of patterns that can explain the entire response, which is interpreted here as 
equivalent to the minimum number of variables that are necessary to construct an index of injury. 

It is worth noting that this two-step hybrid procedure to reduce the problem dimensionality has 
problem-specific advantages over conventional dimension reduction methods such as PCA or Linear 
Discriminant Analysis (LDA). PCA is the optimum method for elimination of collinearities present in 
the data, but is not optimized for class separability, hence will be insensitive to the presence of patterns 
due to elevated injury. LDA is an elegant and "context sensitive" solution, but assumes linearity, 
which may or may not be valid in burn injury. The method employed here incorporates a sophisticated 
model selection process, including nonlinear classifiers such as Artificial Neural Networks. 

Pattern Analysis. The analysis of clustering provides several layers of information. The independent 
pattern analysis performed above provides the ideal number of variables that are necessary to capture 
the injury response. Detailed analysis of cluster membership may reveal further physiological 
information; for instance, in this work we used changing cluster membership between the liver inlets 
(hepatic artery and portal vein) and outlet (vena cava) for any particular metabolite as an indicator of 
altered liver function due to burn injury. 

Variable Elimination and Model Selection. For each model (classifier) type considered, the variable 
elimination task involves finding an optimum list of variables that maximizes the model's prediction 
accuracy. N-fold cross-validation accuracies were used to estimate the real-world accuracy for the 
trained regression models [35,36]: data are separated into N subsets; model training and validation is 
performed N-times, and in each run one of the subsets is used only for validation and the rest for 
model training. N-fold cross validation significantly reduces variability in accuracy estimates due to 
uneven validation data selections as all data is reused for testing the model accuracy. Multiple repetitions of 
the n-fold cross validation ensure that the selection process for the subsets does not affect the results. 

For a given set of variables, comparison of alternative regression models is a straightforward task 
that can be based on the cross-validation results. If the number of variables is preset, the variable 
elimination is also a straightforward, albeit computationally intensive procedure: this problem is a 
combinatorial-optimization problem, where a set number of variables are chosen to maximize cross 
validation accuracy. In this work, we employed the variable selection algorithm in WEKA which 
employs a genetic algorithm to select the best variables to maximize the cross-validation accuracy for 
each classifier. This was followed by a ranking subroutine (Best-first algorithm in WEKA, a greedy 
step-climbing algorithm augmented with backtracking facility) which was used to select only the 
preset number of variables. 

However, the number of variables is typically a confounding factor. Briefly, the training accuracy 
improves in regression as the number of variables/free parameters increases, but this results in 
overtraining, such that the model is not valid for real data {i.e., the model is not generalizable), hence 
cross-validation accuracy will decrease after a certain point. The straightforward approach is to 
evaluate the cross-validation accuracy as a function of the number of variables used in the model, and 
choose the minimum number of variables that provide a high cross-validated accuracy, but this process 
is a very computationally intensive task since the variable elimination described above has to be 
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repeated for each set number of variables, and this process has to be repeated for all possible models. 
Further, there are often multiple good solutions that have very similar accuracies. Therefore, in 
choosing the best combination of model, variables, and accuracy, determination of the proper weights 
of these competing objectives becomes a subjective decision. 

By comparison, the number of independent patterns identified as described above, provides a 
simple criterion that can be determined through an objective and quick process. This approach also 
avoids issues during selection of the best regression model {i.e., the decision of which model to use for 
the index of hypermetabolism), since number of variables employed by each model becomes an a 
priori set quantity that is equal for all tested alternatives. 

3.2. Methods 

Data Preprocessing. The data per rat include metabolite concentrations of each of the PV, HA, and 
hepatic veins, and metabolic fluxes across the liver on day 4 post-burn. To eliminate outliers, variables 
with values beyond the median ± 2 x interquartile range for the group were considered as missing [37]. 
Since individual rat data are used in the construction of a hypermetabolism index, rats with >30% 
missing extracellular metabolite measurements were removed from the experimental dataset. Overall, 
7 out of 12, 7/12, and 7/13 rat datasets were retained for the sham, 20% and 40% TBS A burn 
conditions, respectively. The missing values were replaced by the median of the measurements of that 
group. This resulted in a data matrix of 165 measurements on each of 21 animals. 

Analysis. To identify major patterns in the dose response to burn, k-means clustering was performed in 
MATLAB. The average of each variable was calculated for each burn group (sham, 20%, 40%). Each 
variable was normalized to [-1 1] interval. The Euclidian distance was used as the clustering 
performance which provided better silhouette values compared to other distance measures (results not 
shown). Each clustering was run with >10 replicates and that the cluster means and their membership 
remained similar was confirmed (results not shown). It was observed that the changes in the cluster 
centroids were negligible between clustering runs, an indicator that the number of clusters was well 
chosen and cluster separation achieved was close to ideal. The centroid of each cluster was identified 
as the dose-response pattern. 

Training, cross-validation, and selection of regression models were performed in WEKA data 
mining software [36]. The analysis was performed on per-rat data. The following classifiers were 
tested: Linear Regression (LR), normalized gaussian Radial Basis Function Network (RBFN), Neural 
Network (NN) (multilayer perceptron), Sequential Minimal Optimization algorithm for training a 
support vector Regression model (SMOR), M5P Decision Tree (M5P-DT), Decision Table (DT), M5 
Rules (M5-R). For each method, variable selection was performed as described above. The accuracy of 
each regression model was then evaluated via five 10-fold cross validations, where all randomizable 
variables in the model, as well as the partitioning of data into 10 folds, were randomized. 

4. Results 

Clustering. Analysis of the mean silhouette results demonstrated that increase from four to six cluster 
increased clusters separation only marginally (2.2%), hence four was chosen as the optimum number 
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(Figure 2). Figure 3 displays the results of clustering along with the centroids of each cluster as a 
function of burn injury degree. Table 1 lists all the variables included in the analysis, their averages for 
each burn group, and the cluster membership for each variable. 

Figure 2. Cluster analysis indicates four distinct patterns present in metabolic profiles: 
(A) Number of clusters vs. sum distance from cluster centroid. (B) Number of clusters vs. 
silhouette value. (C) Silhouette values for each metabolite within the four clusters 
identified; a value of 1 indicates the metabolite profile is identical to the cluster centroid; 
-1 indicates the metabolite displays a perfect inverse profile. 




Number of Clusters Number of clusters Silhouette Value 

Figure 3, Identified Clusters. (A-D) Profile of metabolites classified within clusters 1-4. 
(E-H) Centroids (patterns) for each cluster. 
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Table 1. Variables included in analysis and cluster membership. 



Variable 
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Cluster 


a nam 

(median ± IQR) 
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0.2 ±0 
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I 1 A 
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Fraction of Methlyated 
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rlA 
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U ± U 
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Low Density Lipoprotein 
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HA 
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HA 
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Table 1. Cont. 



Variable 
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Total cholesterol/HDL ratio 




i 
i 


Z.UJ ± u.z 


7 14-01 
Z. 1 ± U. 1 


7 7 4-0 
Z.Z ± U 


Valine 




i 
i 


171 £3 4- S4 408 
1Z1.0J ± J4.4Uo 


704 83 4- 74 ft03 
ZU4.0J ± Z4.0UJ 


70ft ft 4- 7 1 QOQ 

zuo.o ± zi.yuy 


A CT/ ATT 
AM/ AL 1 


I I A 

rlA 


7 

z 


7 4- O 7S7 
Z ± U. /jZ 


7/184-1 O/l 
Z.40 ± 1 .U4 


7 07 4- O 733 
Z.U / ± U. / J J 


unionae (^ij 


I 1 A 

rlA 


7 

z 


1 no j. m ? 
luy ± IU.j 


1 1 A 4- 1 Q 

1 14 ± iy 


1 OQ 4- 1 O S 

iuy ± iu.j 


Fraction of Oxyhemoglobin 

(r(Jz.tlu ) 


TJ A 


7 

z 


Q1 9 4- 4 8^ 

y i .z ± 4.o j 


Q1 c + 1 It 
yl.J ± Z. 1 J 


84 8 4- 8 3^ 
o4.o ± o.JJ 


OOUlUm ^INaJ 


UA 
I 1 , \ 


7 


143 4- ft 


1 S3 4- 1 f. s 
1 J J =t 1 O.J 


1 49 4- Q 
1 4Z =t y 


Partial pressure of C02 (pC02) 


HA 


2 


43.95 ±4.9 


51.3 ± 12.65 


46.8 ± 1.4 


pn 


HA 


7 

Z 


7 77 4_ n oso 
/.z / ± u.ujy 


7 3 4- 0 034 
/.J ± U.UJ4 


7 77 4- O 047 
/.Z / ± U.U4Z 


Partial Oxygen pressure (P02) 


HA 
rlA 


7 

Z 


1 Oft 7 4- 1 Q 8S 

iuo.z ± ly.sj 


I 1 7 3 4- 1 8 S 

I I /.J ± lo.J 


1 OO 3S 4- 3S S7S 
1UU.JJ ± JJ.JZJ 


Tyrosine 


I 1 A 


7 


70 814-17 QftQ 


88 1 7 4- ft 78S 
88. 1 / ± O. /OJ 


78 ft/1 4- 8 37S 
/0.04 ± o.jZj 


Acetoacetate 


PV 


2 


115.35 ± 115.91 


134.15 ±91.278 


50.81 ±115.089 


Alanine 


rv 


7 

Z 


3 o 1 £ o _i_ in 77/1 
-591.5s ± 10. /z4 


/I £7 O _l_ 113 /18 

4jZ.y ± 1 1 j.4o 


3HC 77 _l_ flC A/1 

395. zZ ± y5.o4 


Alanine Amino Transferase 

/at r r\ 

(ALL) 


r V 


7 

z 


/I7 4- 1 1 

4Z ±11 


A ft 4- /I S 

40 ± 4.J 


/I7 _i_ /I c 
4Z ± 4.J 


p-hydroxybutyrate 


PV 


2 


81.41 ± 87.551 


103.82 ±64.453 


88.47 ±49.107 


Blood Urea JNitrogen (rsUJNj 


DA / 

r V 


7 

Z 


nice 
1Z ± j.j 


1 A 4- 7 S 

10 ± Z.J 


1 A 4_ 3 S 

14 ± j.j 


Fraction of Oxyhemoglobin 
yrKJztio ) 


DA/ 
r V 


7 

z 


73 ft 4- 1 7 4 
1 J.O ± 1Z.4 


87 7 4- 1 8 7S 
oZ. / ± 18. / J 


73 7 4- 1 S 7 
/J. / ± IJ.Z 


Glucose 


PV 


2 


143 ± 126.5 


183 ±55.5 


147 ±48.5 


Potassium (K) 


DA/ 

r V 


z 


'2 7 4-1 7 S 
J. / ± l.Zj 


/l 1 4_ O 7 

4.1 ± U. / 


3 7 4_ O A 
J. / ± U.o 


Partial pressure of C02 
(puuzj 


DA/ 

r V 


Z 


SS A 4- 3 OS 

JJ.O ± j.Uj 


AO SS 4- S S 

oU.JJ ± J.J 


S7 S 4_ 7 3 
J /.J ± Z.J 


pH 


DA/ 

rv 


z 


7 78 4- n 01 
/.Zo ± U.Ul 


7 3 4- 0 077 
/.J ± U.U / / 


7 7A 4_ O 078 
/.Zo ± U.UZo 


Acetoacetate 


SHVC 


2 


126.21 ±27.78 


126.21 ± 13.89 


30.96 ±95.246 


Aspartate Amino Transferase 
(AST) 


SHVC 


2 


84 ± 16 


85 ± 13 


80 ±31 


Fraction of Oxyhemoglobin 
(F02Hb) 


SHVC 


2 


23.3 ±5.55 


28.7 ± 17.65 


14.1 ± 14 


Leucine 


SHVC 


2 


218.34 ±52.758 


261.92 ±91.466 


210.42 ±147.222 



Metabolites 2012, 2 



468 



Table 1. Cont. 



Variable 


Vessel 


Cluster 


a nam 

( median ± i i< ) 


ZUvc 

(, median ± i k ) 


4Uvc 

(median ± i i^k ) 


pH 


cm/r 


7 

Z 


/.zy ± u.uz 


7 17 4- n oss 

/.JZ ± U.UJJ 


7 7S 4- n 07 1 
/.Zj ± U.UZ1 


Partial Oxygen pressure (P02) 


onVL 


z 


77 7 4_ n 77S 

z /.z ± U. /Zj 


A 1 Q 4- 1 S A 

4i. y ± i j.o 


1 S 1 4- 1 1 A 
1 J.l ±11.4 


Alkaline phosphatase (ALP) 


HA 
HA 


1 
J 


1 07 -1- OA S 

iyz ± yo.j 


111 4-71 
1 1 J ± Z 1 


1 04 4- 74 
1 y4 ± /4 


Aspartate Amino Transferase 
(AL1) 


I 1 A 


1 

J 


AO S 4- 1 4 7S 

4y.j ± 14. zj 


1A 4- 7 S 
JO ± /.J 


10 4- 1 1 S 

jy ± 1 1 . j 


p-hydroxybutyrate 


I l A 
1 1 A 


1 
J 


1 8? £1 4-1 04 804 
ISj.01 ± iy4.Sy4 


AS 4S 4- 7Q 1 17 

0j.4j ± jy. i jz 


1 AS 7 4- 1 04 1S1 
10J.Z ± 1U4.JJJ 


Blood Urea Nitrogen (BUN) 


HA 


3 


12.5 ±7 


13 ± 1.5 


17 ±4.5 


ualcium (uaj 


riA 


i 
j 


8^4-1 7^ 

8. J ± l.Zj 


7 4- 1 

/ ± 1 


8 8 4_ O 8^ 
0.0 ± U.OJ 


Cholesterol 


I 1 A 

HA 


J 


CA C 4_ 1 /I 7C 

jo. J ± 14. / j 


CO 1 1 o 

58 ± 1 J 


71 4_ 7 1 

15 ± 21 


Fraction of Deoxyhemoglobin 
in Total Hemoglobin (FHHb) 


I 1 A 

1 1 A 


1 

J 


n i 4- a a s 

U.J ± 4.4J 


n 7 4- n /i 

U.Z ± U.4 


O 1 4- 7 1 S 

U.J ± /.lj 


Glutamine 


I 1 A 


1 
J 


AHA 4S 4- 07 740 
4U0.4J ± y i .Z4y 


7AA 1 1 4- 47 07 1 
J44.J1 ±4/.UZl 


41 4 14 4- 74 774 
414. J4 ± Z4. / /4 


High Density Lipoprotein 

(I IUL J 


I I A 


1 


JJ 4= J 


77 4-18 
Z / =t 1 o 


1A 4_ A 7S 
JO =t O. / J 


nisLiuine 


I [ A 
I 1 . \ 


1 
J 


as 1 4- 0 404 
OJ.J 4= U.4U4 


A1 07 4_ Q 4S1 

o i .uj =t y.4j J 


A7 1Q 4_ 8 07 

o / . jy 4= o.u / 


.Lactate 


I 1 A 


1 


n si + n Q44 
u.04 =t u.y44 


n ss 4- n Q47 
u.jo ± u.y4z 


1 47 4- 1 40 
J.4 / ± 1 ,4y 


Lysine 


HA 


3 


220.96 ±66.163 


181.93 ±25.467 


224.74 ± 13.38 


Phenylalanine 


I 1 A 

1 1 A 


1 

J 


as no 4- 7 sa/i 
oj.uy ± /.jo4 


A1 7 4-1 SQ8 
OJ.Z ± J.jyo 


A7 17 4-1 707 
O / . 1 / ± 1 .ZUZ 


Total C02 


HA 


3 


21 ± 1.5 


19 ± 1.5 


25 ±3 


l otal rrotem ( l r ) 


riA 


i 
j 


J. J ± U.J 


1 7 4_ O 7 
J. / ± U. / 


A H 4- A 1 

4. / ± U.J 


Alkaline phosphatase (ALP) 


rv 


.5 


1 7/i 4_ ah 


1 ci ill c 

15 1 ± jj.5 


1 Q A 4_ 70 

184 ± /8 


unionde (uij 


DA J 

r V 


i 
j 


11A I 1/1 r 

1 lu ± 14. j 


1 no 4_ 1 1 

iuy ± i j 


117 4-7 

1 1Z ± / 


Fraction of Deoxyhemoglobin 
in i oiai riemogiODin (i^ririD ) 


DA/ 
rv 


1 

J 


10 7 4- 1 1 7S 

ty.z ± ij.zj 


0 7 4-10 7 

y.z ± iy. / 


1 a a 4- 1 s as 

18.0 ± 1 J.8J 


High Density Lipoprotein 
(HDL) 


PV 


3 


36 ± 12 


33 ±4.5 


35.5 ±2.5 


Histidine 


rv 


1 

J 


1 A 0 Q/l 4- 1 A 77 A 
14U.y4 ± lO.ZZO 


170 01 4- Q 87 1 

izu.yj ± y.ozi 


1/17 78 4-10 810 

14/. /s ± iy.8jy 


Lactate 


PV 


3 


1.07 ±0.095 


0.68 ±0.632 


4.38 ± 1.986 


Sodium (Na) 


DA J 

rv 


i 
j 


1 A A 4- A S 
140 ± O.J 


1 A A 4_ 17 ^ 

14o ± IZ.j 


1 A O 4_ 7 S 

i4y ± z.j 


Ornithine 


DA 7 

rv 


5 


i 1 n q 1 4_/i 7^1 
1 19.81 ± 4. /51 


1 f\H 1 A 4_ A 1 C A 

lU/.lo ± 6.154 


1 Tfl 07 4_ 1 7 OA 

izu.sz ± i /.yo 


Proline 


DA/ 

rv 


i 
j 


1A1 774- 10 1^1 

lot .ZZ ± 18. 1 j 1 


1 A7 OA 4_ 8 O 1 1 

loZ.Uo ± o.U 1 j 


171 7 4- IS 1 77 
1 1 Y.Z ± JJ.1Z/ 


Fraction of Deoxyhemoglobin 
in Total Hemoglobin (FHHb) 


on V 


1 

J 


71 1 4- A 
/ 1 . 1 ±0 


SO 4- 1 S 8 

jy ± ij.o 


ao 4- 1 1 4s 

8U ± 1J.4J 


riign uensuy lipoprotein 
(HDL) 


SHVC 


3 


26.5 ± 6 


24 ±5 


33 ±0 


Histidine 


SHVC 


3 


53.01 ±9.45 


52.2 ± 18.937 


61.6 ± 6.314 


Lactate 


SHVC 


3 


1±0.34 


0.74 ±0.141 


3.81 ±2.135 


Lysine 


SHVC 


3 


205.05 ±79.978 


197.38 ± 80.877 


210.04 ±29.84 


Sodium (Na) 


SHVC 


3 


148 ±7 


148 ± 12.5 


149 ±7.5 


Partial pressure of C02 (pC02) 


SHVC 


3 


55.2 ±6.9 


51.3 ±3.75 


67.4 ±7 
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Table 1. Cont. 



Variable 


Vessel 


Cluster 


Cham 

snam 

(median ± IQR) 


ZU Tip 

(median ± IQR) 


4U Tip 

(median ± IQR) 


i oiai luz 


Olivia 


i 


99 + A 

ZZ ± 4 


99 + 9 

ZZ ± z 


9S + H 7S 

Zj ± U. / J 


Tyrosine 


onvt 


i 

j 


A1 88 4- A A8 
Ol.OO ± O.OO 


A 1 90 + 0 771 

oi.zy ± v.i 15 


A8 17 + 71 
OO.J / ± /.J 


Acetoacetate 




A 
4 


1 1 A 1Q 4-11 8 1 Q 
1 10. Zy ± ZZ.oiy 


1 1 n 14 + 47 A91 
1 1U.J4 ± 4 /.OZj 


70 AS + 1 1 740 
/U.OJ ± J 1 . /4y 


Alanine 


HA 


A 

4 


1/1 Q 11 4- A 1 Q^/l 


9Q7 19 + 99 71 S 
Ay /.JZ ± ZZ. / 1 j 


981 A + 11 081 
Zo 1 .4 ± J J.UOJ 


Albumin 


HA 


A 

4 


1 .yj ± U. 1 / J 


i 9 + n i 

1 .Z ± U.J 


114-0 9S 
1 . J ± U.Zj 


Alanine Amino Transferase 




A 
4 


AQ ^4- 14 OS 
Hy.J =t 14. Zj 


1A + 7 S 
JO ± / . J 


1Q 4- 1 1 S 
J7 ± 1 1 . J 


/\IIlIIlUIlld 


HA 


A 
H 


AQ 7 4- A ^11 


17 1 + 4 081 
J / . J + 4.Uo 1 


91 1 + 4 A1 A 
Z J .1 + H-.O 1 o 


Asparagme 


HA 


A 
4 


41 Q8 4- ^ 488 
4 1 .70 ± j.4oo 


14 S 1 +1 Q4S 
J4.J 1 ± J.y4j 


19 71 4- A S8 
JZ. / J ± o.Jo 




HA 


A 
H 


11 81 4- 1 S4Q 
1 1 .OJ =t 1 . J4y 


Q 9S + S 8S9 
y.Zj =t j.OJZ 


10 1 A + 9 A17 
10. 10 + Z.OJ / 


riemogiODin \/o) 


HA 
1 LA 


A 
4 


19 4+1 I s ; 

1Z.4 =t 1 Jj 


19 + 94 
1Z =t z.4 


1 0 4 + 1 OS 
1 U.4 + 1 .Oj 


r 1UW IvdlC 


HA 
1 LA 


A 
H 


n as + n 1 4^ 

U.OJ =t U. 140 


n s + n nAQ 

O.J + O.OO" 


0 S + 0 49A 
O.J + 0.4ZO 




HA 
[ la 


A 
H 


1 AQ + 87 
1 DV 4i 0 / 


1 14 + 4A 

1 JH- 41 H-O 


1 4A S + 1 A 

1 H-O. J + 1 O 


uiuiamaie 


HA 
HA 


A 
4 


7A 8S + A AA9 
/0.OJ ± O.OOZ 


S4 1 S + 1 Q 1 17 
J4. lj± ly.lj/ 


4A 84 + 1 7 4S1 
40.04 + 1 /.4JJ 




HA 


A 
H 


9 1 A 7S + 1 A 1 OS 
ZIO./J 4= lO.lyj 


1 84 88 + 9S Q74 
lo4.oo + Zj.y/4 


1 A1 7A + 1 Q 991 
1 OJ. / 0 + 1 7.ZZJ 


ncITlalOClTl 


HA 


A 
4 




IS + 7 


11+1 

J 1 + J 




HA 


A 
H 


ill 89 +47 491 


1 1 9 OS + 1 Q 1QA 

i iz.uj + iy.jyo 


1 08 01 + 19 S9Q 
100.01 + jA.jAy 


U/a+oooi i lin / It 1 

roiassium (is.) 


HA 
HA 


A 
4 


4 9^4-9 1 s 
4.ZJ ± Z.lJ 


4 9+1 4S 
4.Z =t 1 .^O 


1 Q 4- 0 7S 
j.y ± U. / J 


leucine 


HA 
1 LA 


A 
4 


1AQ 7 + 1 74 8 1 7 
joy. / =tl/4.ol/ 


97S 1 + 19 QSA 
Z/j.l =t JZ.yjO 


98S 9 + 87 QQQ 
ZOj.Z + o 1 .yyy 


Methionine 


HA 


4 


52.33 ± 11.417 


40.59 ± 8.321 


38.75 ±9.047 


Proline 


H A 


A 

4 


1 A/l S9 + Q OS A 

io4. jz ± y.yjo 


1/17 74- 97 A Q9 

14 /. / ± z /.4yz 


1/17 AQ + 7 111 
14 / .Oy ± /.JJJ 


serine 


HA 
I 1 A 


A 

4 


991 81 + S4 1 OA 
ZZj.OJ ± j4.iyo 


1 SO SI + 1 9 879 

ijy.jj ± iz.o /z 


1 AQ OA + 9A 87 

loy.uo ± ZO.O / 


i ngrycenaes ^ioj 


HA 


A 

4 


(.1 4_ IS 7S 
OZ ± J J. / J 


90 + 7 

zy ± / 


9A + A 
ZO ± 0 


Threonine 


HA 
1 1 A 


A 

4 


9SA HI + 1A 78S 

Zjo.Uj ± jo. /oj 


990 07 + 41 119 
zzy.u / ± 4 1 . J JZ 


909 04 + 1 S 004 
ZUZ.U4 ± 1J.UV4 


Valine 


HA 
I 1 A 


A 

4 


914 08 + 1 07 898 
Zj4.yo ± 1U /.oZo 


1 84 IS + 1A SSQ 
1 o4.J J ± JO.J jy 


1 78 74 + 11 8S 1 
1 /o. /4 ± JJ.OJ 1 


Very low density lipoprotein 

A/T FIT 4 


HA 
1 LA 


A 
H 


1 1 4_ n s 

1 J 4= U. J 


A + 1 
0+1 


S + 0 S 

J + O.J 


A 1 i i in 1 n 

/\iuuniin 


PV 
r V 


A 
4 


9 4-0 97S 
Z ± U.Z / J 


1 i 4_ n is 

1 . J ± U. J J 


114-09 
1 . J ± U.Z 


/\IIlIIlUIlld 


PV 
r V 


A 
H 


1 f)8 AS + A1 S9S 
lUo.OJ =t 01. jZj 


41 74 + 11 191 
4J. /4 + JJ.1Z1 


79 A7 + 1 0 70Q 
/Z.O / + 1 U. / Vy 


A or^o t*o m y\ c± 

/\&pdrdgine 


PV 
r V 


A 
H 


S 1 f)8 + A 88A 
J 1 .UO =t O.ooO 


19 OQ + 1 80S 

JZ.U7 4= J.07J 


11 74 + 0 817 
J 1 . /4 + U.OJ / 


A ct^ortotp 

/\&pdridic 


PV 
r V 


A 
H 


10 18+11 f)Q7 
1 y . 1 o =t ii .uy / 


19 7 4-1 SQ4 
1Z. / + i .jyM- 


19 OA + 9 988 
1Z.UO + z.zoo 


i^diLiurn v^dj 


PV 
r V 


A 
H 


8 7 + 9 9S 
O. / =t Z.ZJ 


8 i + n ss 

O.J + U.JJ 


70 + 1 AS 
/ .y + 1 .0 J 


riemogiODin ^ /oj 


PV 
C V 


A 
4 


1 4 A + fl QS 

14.0 ± u.yj 


114-1 7S 
1 J ± 1 . / J 


19 8 + 0 S 
IZ.O + u.J 


l T ii torn o 


PV 


A 
H 


A1 1 A + 11 978 
OJ. 10 4= 1 J.Z / o 


49 99 4- 1 1 1 48 

H-Z.ZZ + 11.14-0 


48 44 + 1 S 1 1 Q 

Ho.HH + 1J.J17 


It 1 1 1 torn ihp 
U 1 U LdllllllC 


PV 
r V 


A 
H 


970 1 ? + SO 01 1 


9si ns + io 1 no 

ZJJ.OJ + /U.1U7 


9S0 ?Q + 33 943 
z*j\j.z*y + jj.zh-j 


Glycine 


PV 


4 


273.49 ±48.285 


222.25 ±33.939 


204.59 ±46.681 


Hematocrit 


PV 


4 


43 ±3.5 


40 ±5 


38 ±2.5 


Methionine 


PV 


4 


44.77 ±4.889 


44.17 ±4.749 


43.82 ±3.209 


Partial Oxygen pressure (P02) 


PV 


4 


69.9 ± 11.75 


54.2 ± 8.1 


62.9 ± 15.65 


Serine 


PV 


4 


202.59 ± 14.805 


148.88 ±30.334 


156.79 ± 14.701 


Triglycerides (TG) 


PV 


4 


52 ±43 


25.5 ± 11.75 


23 ±2 
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Variable 


Vessel 


Cluster 


CU ftl « 

dnam 

(meaian ± ii^ik) 


( iviArliitn _l_ ¥ /" \D \ 

(mcciian ± ii^k ) 


4Uvc 

(meai an ± ii^K) 


Total Protein (TP) 


PV 


4 


5 ± 1.25 


4.7 ±0.9 


4.8 ±0.85 


Very low density lipoprotein 
f\n fit 


T>\7 
r V 


A 
4 


1 1 4- s 


f. 4- 1 
0 ± 1 


J ± U.J 


Alanine 


onvt 


A 

4 


7/1A 18 4- SI 8Q1 
Z40.J0 ± JJ.oyi 


1 S A 77 4- A 7 S7Q 
1 JO. /Z ± 4Z.jZy 


1 S7 7^4-1/1 07 
IjZ.ZO ± 14. UZ 


Albumin 


orlVU 


A 

4 


I .y ± U.ZJ 


1 1 4_ 0 7S 
1 .J ± U.ZJ 


17 4- 0 7S 
1 .Z ± U.ZJ 


Alkaline phosphatase (ALP) 




A 
4 


1 AA -I- A7 S 
100 ± 0 /.J 


1 7Q 4_ 7 1 

izy ± zi 


1 47 4- 87 
14 / ± OZ 


Alanine Amino Transferase 


Oil V L 


A 
H 


4fl -1- 1 S S 
4U =t 1 J.J 


jo =t 14 


37 4- S S 

J / 4= J. J 


/\mmoma 


oil V 


A 
4 


40 AA 4- 7 A77 
4U.44 =t Z.OZZ 


71 86 4- 1 0 QS1 

zi.so ± iu.yj i 


1 S 01 4- Q 87S 
1 J.U1 ± y.o / J 


Asparagine 


SHVC 


4 


27.87 ± 1.556 


20.1 ±3.165 


18.34 ±0.097 


Aspartate 


oriVU 


A 

4 


1 1 57 4_ 1 Q1Q 

i j.j / ± i.yjo 


1 0 17 4- /I Q17 

iu.j / ± 4.yjz 


8 /174- 0 817 
8.4 / ± U.O 1 / 


p-hydroxybutyrate 


orlVU 


A 

4 


77S AO -1- 4S7 1 1 1 

z / j.oy ± 4j /.j 1 1 


7^7 78 4- 744 001 
ZJ /.Zo ± Z44.UU1 


771 SI 4- 1 74 Q44 
ZZj.jl ± 1 /4.V44 


uaicium (\,& ) 


crn/ri 


A 

4 


8 A 4- 7 S 
O.O ± Z. J 


7 O 4- O 1 
/ .y ± U.J 


8 4- 1 S 
o ± 1 . J 


unionae (i^ij 


oil VL 


A 

4 


111 4-1/1 s 
1 1 1 ± 14. J 


1 1 n s 4- a i< 
1 IU.J ± 4. / J 


1 OQ 4- 1 1 

iuy ±11 


Hemoglobin (%) 


orivi^ 


4 


1 7 7S 4_ 0 87S 
1Z. / J ± U.o / J 


1 O 7 4- 1 7S 
1U. / ± 1. /J 


11 A 4- 0 AS 
1 1 .0 ± U.Oj 


Fraction of 

1 of nr\v\;nPtvir\(r mnin / h 1 1 1 1— lr\i 

i^drDUAyncinugiuuiii \ r laji idj 


Oil V L 


H 


7Qf) 4- 73 7S 

Z;?U =t / J. / J 


1 30 4- 94 
1 JO 4i ZH- 


1 lf\ 4- SO 
1 / O =t JU 




Oil V ^ 


A 
H 


963 4- 73 
ZOJ =t / J 


1 8^ 4- 79 S 
1 oO =t / Z.J 


1 S9 4- 3A 9S 
1 JZ =t JO.Z J 


UlUldlllalC 


oil V L 


A 
H 


If. 77 4_ AT. Ill 
/O. / / =t 4J.Z / / 


4Q fj~) 4_ 1 S 7Q 
47.0Z =t Lj./y 


44 7 4- 1 Q 0Q3 


oiuiairiine 


oil v 


A 
4 


7QQ 7 + « 737 
Lyy. I ± JO. / J / 


7S4 87 4- Q7 78A 
ZJ4.5Z ± y / . / 00 


778 OA 4- Q3 784 
Z.Z.O.V0 41 yj.Zo4 


urycme 


oil v 


A 
4 


1 Q7 AS 4- S7 864 
iyZ.4j ± j /.004 


1 48 04 4- 8 773 
140. U4 XO.//J 
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Phenylalanine 


SHVC 


4 


54.68 ±3.359 


52.83 ± 15.86 


49.57 ± 14.92 


Pro I mp 


Oil V v 


4 


1 37 87 4- 37 486 


1 1 < ^4 _|_ f. 039 


113 16+16 448 


Serine 


SHVC 


4 


166.82 ±20.803 


105.99 ± 12.727 


92.92 ± 42.759 


Triglycerides (TG) 


SHVC 


4 


25 ±28.5 


20 ±3 


20 ±0 


Threonine 


SHVC 


4 


198.18 ±20.372 


186.17 ±59.461 


157.26 ±6.575 


Total Protein (TP) 


SHVC 


4 


4.9 ± 1.05 


4.3 ±0.45 


4 ±0.85 


Very low density lipoprotein 
(VLDL) 


SHVC 


4 


9.5 ±4.5 


5±0 


5±0 



Identification of Major Patterns. While there are four major patterns in the metabolic response to burn, 
they are not necessarily independent. For instance, clusters 1 and 4 apparently display the inverse of 
the same response. To identify the number of independent patterns, SVD Analysis was employed. This 
analysis revealed that two patterns explained 98.6% of the total variability observed, indicating that the 
remaining two patterns were explainable within a margin of error of <2%. The patterns highest 
weighted in SVD were clusters #2 and #4. It should be noted that the pattern selection here is based 
purely on overall trends, unlike the analysis of Yang et al. [34] which identifies critical time points of 
sudden changes in gene expression. Since the resolution of data for burn injury is limited (i.e., we have 
only three levels of burn injury), it is not possible to employ a similar method to identify the degree of 
injury where the post-burn hypermetabolic response kicks in. 
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Analysis of Cluster Membership. In general there were no strong trends observed in the variables 
comprising the clusters (Table 2). Cluster 2 had nearly equal membership from all vessels. Cluster 4 
had the most variables (68 total). Albumin, asparagine, aspartate, glutamate, methionine, serine, 
hematocrit, triglycerides and VLDL concentrations systemically (i.e., in all vessels observed) were 
selected to cluster 4, which displayed a general decreasing response with increasing burn. pH and 
FC^Hb were in cluster 2, which showed a peak at 20% burn. Lactate had increased systemically and 
was in cluster 3 with a sharp increase at 40% burn; HDL and histidine also belonged in this cluster 
systematically, although the sharp increase at maximum burn was not present. LDL and cysteine 
increased systemically, and were in cluster 3. 



Table 2. Cluster membership as a function of vessel/flux. 





# variables in vessel/flux 




Hepatic 


Portal 


Vena 


Liver 


Total 


Cluster 


Artery 


Vein 


Cava 


Flux 




1 


9 


20 


15 


9 


45 


2 


8 


10 


6 


8 


26 


3 


15 


9 


9 


15 


36 


4 


23 


16 


25 


23 


68 


Total 


55 


55 


55 


5 





To analyze the role of the liver in the dose response to burn injury, the variables that were selected 
to the same cluster in the liver inlet (PV and HA) but different outlet (SHVC) were identified. FCOHb 
(HA and PV in cluster 1, SHVC in cluster 4), pC02 (2/3), ALP (3/4) were the three variables identified. 

4.1. Model Selection and Index of Burn Injury Severity 



The variable elimination process was performed with each model with the number of variables set 
at two. The two best models were the M5-Rules and NN, with very similar cross validation accuracies 
(Table 3). The formula for the index of burn injury severity developed via the NN model is: 



Burn(%mSA) = 2.11 



1 



-(-3.08 VLDL(HA)-4.63 AcetoAcetate (SHVC)-l .55) 



yl + e~ 

Table 3. Comparison of tested models. 



-1.0388 





Average Relative Accuracy 


Variable 1 


Variable 2 


RBF 


46.06 ±4.59 


Ammonia (HA) 


Ornithine (SHVC) 


DT 


70.54 ±7.83 


Asparagine (SHVC) 


ALB (HA) 


M5P-DT 


72.36 ± 1.03 


ALB (HA) 


Ammonia (HA) 


LR 


73.28 ±0.55 


VLDL (HA) 


Acetoacetate (SHVC) 


SMOR 


75.92 ±2.20 


VLDL (HA) 


Acetoacetate (SHVC) 


M5-R 


87.82 ± 1.31 


VLDL (HA) 


Asparagine (SHVC) 


NN 


87.89 ±2.79 


VLDL (HA) 


Acetoacetate (SHVC) 



An interesting observation was that very low density lipoprotein (VLDL) (cluster #4) was a 
common choice in all the high-accuracy models. The acetoacetate level in the SHVC was the second 
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selection in the NN model, whereas SHVC asparagine concentration was selected in the M5P-R model. 
The top two regression models, NN and M5P-R, had nearly identical -88% cross validation accuracy, 
which has value as a diagnostic index, but with significant room for improvement. 

The effects of using additional variables in the index were evaluated for NN and M5P-R (Table 4). 
The NN model was selected as the regression model of choice since accuracy increased significantly 
(23%) with up to 4 variables. By contrast, the accuracy of the M5P-R model was slightly increased 
with a third variable, but later decreased with addition of variables (it should be noted this is indeed the 
expected situation with proper cross-validation, which can account for over training; without 
validation data the accuracy would have increased monotonously with addition of new variables). The 
predictions of the NN model on a case-by-case basis are displayed in Table 5. However, note that these 
are the results of training on the full set of data, hence the overall accuracy is significantly higher than it 
would be expected in an actual application, which is what the cross validation results in Table 3 report. 



Table 4. Effect of number of variables on regression accuracy. 



Method 


Measure 


2 Variables 


3 Variables 


4 Variables 


5 Variables 




Average Relative 
Accuracy 


87.89 ±2.789 


89.46 ± 1.379 


90.69 ± 1.508 


89.18 ± 1.552 


Selected 
Variables 


VLDL (HA), 
Acetoacetate 
(SHVC) 


VLDL (HA), 
Acetoacetate 
(SHVC), 
p02 (SHVC) 


VLDL 

(SHVC), 
Acetoacetate 
(SHVC), 
p02 (SHVC), 
C02 (HA) 


VLDL (SHVC), 

Acetoacetate 

(SHVC), 

p02 (SHVC), 

Asparagine 

(SHVC), 

C02 (HA) 


M5-R 


Average Relative 
Accuracy 


87.82 ± 1.312 


88.71 ±0.983 


87.66 ± 1.398 


87.11 ± 1.038 


Selected 
Variables 


VLDL (HA), 

Asparagine 

(SHVC) 


Albumin (HA), 
Asparagine 
(SHVC), 
BUN (PV) 


VLDL (HA), 
Acetoacetate 
(SHVC), 
Albumin (HA), 
pC0 2 (SHVC) 


VLDL (HA), 

Acetoacetate 

(SHVC), 

Albumin (HA), 

pC0 2 (SHVC) 

Asparagine 

(SHVC) 



Table 5. Predictions of neural network model (Training on all data). 





Burn % 
(actual) 


Predictions 


Instance 


2 Variable 


3 Variable 


4 Variable 


5 Variable 




Model (%burn) 


Model (%burn) 


Model (%burn) 


Model (%burn) 


1 


0 


0.132 


0.324 


1.35 


1.271 


2 


0 


-0.281 


-0.373 


-0.423 


-0.320 


3 


0 


3.003 


0.109 


0.05 


-0.195 


4 


0 


-0.614 


0.415 


-0.632 


-0.52 0a 


5 


0 


0.309 


-0.713 


-0.011 


-0.013 


6 


0 


-0.416 


-0.213 


-0.451 


-0.635 
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Table 5. Cont. 





Burn % 
(actual) 


Predictions 


lIlSldllLC 


z v aria Die 


j v aria Die 


t v aria Die 


3 v aria Die 
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40 870 

4U.O / U 


17 


40 


40.709 


40.091 


40.045 


39.633 


18 


40 


41.238 


39.938 


39.866 


39.709 


19 


40 


41.221 


38.822 


39.612 


39.620 


20 


40 


39.865 


39.757 


38.984 


38.333 


21 


40 


36.667 


39.739 


39.708 


39.546 


Relative Absolute 


8.49% 


2.54% 


4.10% 


5.04% 


Error 













Finally, since selection of acetoacetate is highly significant as an indicator of mitochondrial redox 
potential, we also tested the prediction success when acetoacetate/p-hydroxybutyrate ratio 
(a commonly used predictor of redox potential) is used [38,39]. Briefly, we repeated the variable 
selection/training algorithm for the following set: acetoacetate/p-hydroxybutyrate ratio (SHVC), 
acetoacetate/p-hydroxybutyrate ratio (HA), VLDL (SHVC), p0 2 (SHVC), C0 2 (HA), as well as 
acetoacetate/p-hydroxybutyrate ratio (SHVC) alone with NN and M5-R models. On its own, the 
acetoacetate/p-hydroxybutyrate ratio was extremely unsuccessful in predicting %TBSA, (average 
relative errors of 110.87 ± 7.444 and 110.26 ± 3.063 with NN and M5-R, respectively). Inclusion of 
other variables marginally improved the results but nowhere near the other results in Table 4: The 
average relative errors were 76.26 ± 15.887 and 64.55 ± 10.283 for NN and M5-R, respectively. 

5. Discussion and Conclusions 

The most important finding of this work is the strong correlation between VLDL levels and the 
extent of burn injury. Prior studies indicate that in burn trauma, VLDL secretion is impaired [40]. This 
is consistent with our observations that for increasing burn, there was a decrease in systemic VLDL. 
Interestingly, VLDL, LDL, HDL and TG are four of the relatively few variables that were systemically 
altered (i.e., belonged in the same cluster independent of the measured blood vessel). However, the 
model developed here suggests that VLDL is the most preferred variable for building a predictor of 
severity of burn injury; ahead of other potential metabolic targets such as hyperglycemia which is 
correlated to insulin resistance [9], alanine/glucose (a systemic cycle for conversion of skeletal muscle 
to glucose during hypermetabolism), or even glutamine and arginine, which tend to be depleted after 
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burn injury [5,41,42]. Of note here is that as clustering analysis shows, other potential variables 
identified in Table 4 demonstrate altering {i.e., belonging to different clusters based on measured 
vessel) injury dose-response profiles; therefore, predictions based on VLDL are least likely to be 
affected by tissue specific variations. Very likely, the fact that VLDL decrease is systemic also renders 
it a more accurate indicator of hypermetabolism. 

The second variable of interest, which was selected commonly for most regression models 
(including M5P-R with 3 or more variables) is acetoacetate (SHVC), which in the context of burn is 
most important as the precursor for P-hydroxybutyrate. While ketone bodies are not a subject of 
intense focus in hypermetabolism, it was previously suggested that the acetoacetate/p-hydroxybutyrate 
ratio reflects the mitochondrial redox potential in the liver [43,44], and in burn patients, a decrease of 
plasma acetoacetate to P-hydroxybutyrate ratio indicates mitochondrial dysfunction and correlates with 
developing multiple organ dysfunction [45]. Our rat data show a significant decrease in acetoacetate in 
the 40% TBSA burn group (Table 1). P-hydroxybutyrate was also decreased, although to a lesser 
extent, such that the acetoacetate to P-hydroxybutyrate ratio was decreased as well. In general, the 
decrease in total ketone bodies is regarded as an indicator of a limitation in oxidative phosphorylation 
after injury [46]. Concomitantly, other variables linked to mitochondrial activity or respiration, such as 
venous CO2, were also selected as additional injury indicators by the Neural Network model. It should 
be noted however, that the acetoacetate/p-hydroxybutyrate ratio by itself was not a particularly useful 
indicator of extent of burn injury. Further, replacing acetoacetate by the ratio severely decreased the 
accuracy achievable, with 100% error range, indicating that since the ratio was unsuccessful both on 
its own and when replacing acetoacetate concentrations, this reduction in accuracy is likely to be at 
least partially due to amplification or measurement noise. Note that P-hydroxybutyrate was one of the 
measurements with highest standard deviations relative to group means. 

A third variable that was commonly selected was asparagine. Asparagine/aspartate is one of the 
cycles that carry amino acids from the muscle tissue to the liver during injury and was previously 
reported to be altered in rat models of burn injury [47]. While asparagine was able to replace 
acetoacetate in the M5P-R model, for the NN model that was ultimately more successful and could 
exceed 90% accuracy acetoacetate proved a better variable to include in the index. 

It is worth noting that VLDL (SHVC) is not significantly different between the 20% and 40% burn 
groups, but is significantly reduced compared to the sham group. It can be argued that the ability to 
differentiate between sham and burn provided by VLDL complements acetoacetate (SHVC), which is 
different for the 40% burn group, but similar for the sham and 20% burn groups. The interesting note 
here is that neither of these variables is simply directly correlated to extent of burn, hence use of 
multiple variables is necessary. This supports our previous findings [21,48] that there are significantly 
different responses observed between the 20% and 40% burn groups, either because the response at 
20% TBSA is significantly less, or possibly as the 20% group is displaying a switch from 
hypermetabolism to a healing/normal phase as early as 4 days after burn injury. This response is likely 
the reason that a nonlinear model, such as the multi-layer perceptron, has higher accuracy than other 
linear models we tested in this study. This result may also justify why hypermetabolism diagnosis is 
still best based on clinical observations because this ad hoc method allows the physician to account for 
the nonlinear behavior based on past experience. 
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Using the 2-variable minimum identified by clustering, an index of burn injury severity was 
developed. The cross-validated accuracies of the best regression model, artificial Neural Network, was 
88%. We also tested the inclusion of additional variables into the regression model. As displayed in 
Table 4, the accuracy of the multilayer perceptron model could be increased to up to 91% with the 
addition of arterial total CO2 (cluster 3) and venous oxygen (cluster 2). While the NN therefore did not 
include any variable from cluster 1, since cluster 4 displays the inverse dose-response to cluster 1 very 
closely, this confirms the previous analysis that 2 clusters are sufficient in capturing most of the 
metabolic response to increased burn injury. 

From an application perspective, VLDL is in the same cluster systemically; hence point of 
measurement is unlikely to affect results significantly (repetition of 2-variable NN model with SHVC 
VLDL use led to only one rat in the sham burn group being significantly misclassified as 20% burn 
animal). Acetoacetate also displays a generally decreasing trend in all vessels. It should be noted that 
for an actual clinical index, ideally all measured variables will be systemically in the same cluster, as 
well as practical to measure. Most variables meet the criteria selected in the regression models. It is 
also likely possible to construct a regression model to predict unsuitable metabolites from more 
easy-to-measure ones, which we have not investigated here; a proper study of such an approach would 
require repetition of the experiments with tail-vein blood sample data complementing the HA, PV and 
SHVC samples to test if these less invasive samples correlate accurately to the data in this work. 

To our knowledge this is the first attempt at creating a quantitative index of burn injury severity; 
however, it is important to realize the limitations in clinical applications. As the animals were not 
subjected to any intervention following burn injury, there were no potential confounding effects from 
nutritional supplementation, which may not be the case in human patients. Obviously, differences 
between rat and human metabolism, as well as effects of age and gender differences on the hypermetabolic 
response [48], will have to be considered for the development of a clinically applicable index. 

These indicators of burn injury may provide insight and clues to new metabolic targets for therapy. 
In addition, the development of a quantitative score to identify the degree of hypermetabolism can 
ultimately provide a practical way to measure the patient response to injury and treatment. As the 
current diagnostic criteria are based purely on clinical observations, patient-to-patient variation in 
metabolism introduces a significant and undesirable degree of uncertainty in the care of the burn 
patient that could be avoided by such a quantitative index. 
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